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(-H ' Abstract 

' A convection problem with temperature-dependent viscosity in an infinite layer is pre- 

I , sented. As described, this problem has important applications in mantle convection. The 

^ ' existence of a stationary bifurcation is proved together with a condition to obtain the critical 

^\ ' parameters at which the bifurcation takes place. For a general dependence of viscosity with 

0^ ■ temperature a numerical strategy for the calculation of the critical bifurcation curves and 

' the most unstable modes has been developed. For a exponential dependence of viscosity 
on temperature the numerical calculations have been done. Comparisons with the clas- 

. sical Rayleigh-Benard problem with constant viscosity indicate that the critical threshold 

^-H ' decreases as the exponential rate parameter increases. 

oo : 

o. ...... 



1 Introduction 

A classical subject in fluid mechanics is the problem of thcrmoconvective instabilities in 
fluid layers driven by a temperature gradient [H 1161 [T]. Gravity and/or capillary forces are 
responsible for the onset of motion when the temperature difference exceeds a certain threshold. 
The existence of a stationary bifurcation for the gravity driven case (the classical Rayleigh- 
Benard problem) has been rigorously proved in [H [TT] . Above the threshold convective solutions 
consist of rolls [TU] ■ Most studies on the Rayleigh-Benard problem consider a constant viscosity 
[131 M US] J although interest on convection problems with temperature dependent viscosity has 
increased [TU IH [171 |6l [71 (TU [20] due to the fact that this dependence is a fundamental feature 
of mantle convection. The formation of rocks indicates that viscosity in the interior of the Earth 
strongly depends on temperature, and this influences the thermal evolution of the mantle [SJIS]. 
Some theoretical or numerical studies with temperature dependent viscosity in flnite boxes can 
be found for instance in Refs. [THIITI [TU [T^] . In [TB] several types of dependence have been used, 
but the most common is an exponential dependence [131 [7] . 

This paper plans a theoretical and numerical study of the instabilities in a gravity driven 
convection problem with viscosity dependent on temperature. The domain, which is finite in the 
vertical coordinate and infinite in the horizontal plane, contains a viscous fluid which is heated 
uniformly from below. In this work, for a general dependence of viscosity on temperature, we 
prove the existence of a stationary bifurcation as the temperature gradient increases. And a 
condition to obtain the critical parameters at which the bifurcation takes place has been found. 
For a viscosity with an exponential dependence on temperature both the critical bifurcation 
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curves and the most unstable modes have been numerically computed. Comparisons with the 
classical Rayleigh-Bcnard problem with constant viscosity indicate that the critical threshold 
decreases as the exponential rate parameter increases. We find the critical bifurcation curves for 
this case as a function of the exponential parameter. 

The paper is organized as follows. In section 2 the formulation of the problem with the 
equations and boundary conditions (be) and the conductive solution are presented. Section 3 
develops the linear stability of the conductive solution and deduces the theoretical results on the 
existence and conditions of bifurcation. In section 4 numerical results on eigenmodes and critical 
curves are presented. Finally in Section 5 conclusions are detailed. 



2 Formulation of the problem 



The considered physical set-up is shown in figure 1. A fluid layer of depth d {z coordinate) is 
placed between two parallel boundless plates. On the bottom plate a temperature To is imposed 
and on the upper plate a temperature Ti = T^^AT = To'~(3d, where /3 is the vertical temperature 
gradient. The governing equations are the continuity equation 



V-v = 0, (1) 



the Navier-Stokes equations 



dtv + {v-\/)v = -—ges - — VP + div{v{T) ■ (Vw + (VC)*)), (2) 

Po Po 



and the heat equation, 



dtT + V-VT = KAT. (3) 



Here v = (vxTVy^Vz) {x,y,z,t) is the velocity vector field, P is the pressure, T the temperature 
field, p is the density, which dependence on temperature is assumed to be p = — a (T — Ti)), 
being pi the mean density at temperature Ti, a is the thermal expansion, K is the thermal 
diffusion, g is the acceleration of the gravity and 63 is the unitary vector in the vertical direction. 
i^(T) is the viscosity that is assumed to be a positive, convex and bounded function of temperature 
T. 

Now, the boundary conditions are introduced. The temperature on the bottom plate, z — 
—d/2, and on the upper plate, z = d/2, are respectively 



T^z^--j^n, T^z.-j.r,. (4) 

The boundary conditions for the velocity correspond to rigid, no-slip at the bottom plate and 
non-deformable and free-slip at the upper plate. They are expressed as follows. 



Vy ^v^ =0 on z = dzVx = dxVy = = on z = ^. (5) 



The simpler solution of the hydrodynamic equations is the conductive solution. For that solution 

the temperature only depends on the vertical component and the fluid is at rest, 

d^T = dyT ^dfT^O, 0. (6) 

Under these assumptions the hydrodynamic equations become, 

= -pgeh-VP, (7) 
= KAT, (8) 
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which solution is the conductive state 



To-[i\z + - 



Pb [z) = Po- Pog 

Vb = 0, 
where Po is an arbitrary constant. 



z 1 



a[5d\ aPz'^ 



(9) 

(10) 
(11) 



3 Linear stability of the conductive solution 

Now, we are interested in the study of the hnear stabihty analysis for the conductive solution. 
Therefore, we perturb it as follows. 



T{x,y,z,t) = Tb{z) +e{x,y,z,t), 
P{x,y,z,t) = Pb{z) + SP{x,y,z,t), 
v{x,y,z,t) = Vb + u{x,y,z,t), 



(12) 
(13) 
(14) 



and introducing these expressions into the hydrodynamic equations the following equations for 
the perturbation fields are obtained. 



V • u = 0, 

dtu + (u • V) It = aege3 - —V (SP) + div(i^(T) • (Vu + (Vw)*)), 

Po 

dt9 + {u-W)9 = KA9 + /3u^ . 

Taking into account the following expressions, 

i^iT) = v{Tb + 9) = u{Tb) + y'{Tb)9 + 0{9'^), 
u{Tb)^v{To- p{z + d/2))^u{z), 
div(zy(T) • (Vu + (Vm)*)) ^Lu + N{u, 9), 

where Lu is the linear part and N{u, 9) includes the nonlinear terms as follows 

Lu = div {v{z) ■ (Vu + (Vw)*)) , 
N{u,9) = A:iY{{v'{Tb)9 + 0{9'')) ■ {Vu+{WuY)) . 

Linearizing eqs. (jlSp the following system of equations is obtained 

V • u = 0, 

dtU = a9ge3 ~ —V (SP) + Lu, 

Po 

dt9 ^ KA9 + f3u,, 



(15) 



(16) 



where displaying the operator Lu = dz {i'{z)dzu) + v{z)/S.' u -\- v\z)^Uz and A' is the laplacian 
operator in the variables x and y. 

The system (|16|) together with its boundary conditions may be expressed in dimcnsionlcss 
form defining new variables = r/rf, t' = Kt/tP, u' — du/K,p' = d'^p/ {poKvo), 9' = {9—Ti)/ [3d. 
After rescaling the variables and dropping the primes the system is left as follows. 



V • u = 

—dtu = 
Pr 

dt9 ^A9 + u 



V (dP) + Lu + ROes, 



(17) 
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Here, the Raylcigh number is defined as i? = d^agP / {vi^K) and the Prandtl number is Pr 
uq/K, where vq = v{—l/2). 

The dimensionless boundary conditions for fields in p7p are 



— u,, 



on z 



-1/2, e ^ d^u^ ^ d^Uy 



0onz = l/2. (18) 



3.1 Decomposition into normal modes 

As the domain is infinite in the plane, the eigenvalue problem can be solved by using the 
normal modes method. For each mode a separate variable solution can be found S{z) ■ /(x, y, t), 
where f{x,y,t) can be expressed as a normal mode f{x,y,t) = ^'^Q')'t+^ik:oX+kyy) _ Therefore we 
look for solutions as follows, 



( Ux \ 



V 5P J 



( Ux{z) \ 

Uy{z) 

W{z) 

e{z) 
V M^) J 



T(fe)t+i(fcxX + fcyJ/) 



(19) 



where cr(fc) are the eigenvalues, that initially can be complex numbers. It represents the growing 
rate of the normal mode at the wave number {kx, ky) ~ k. If the solution belongs to in kx^ ky, 
the inverse Fourier transform of the normal modes can be considered the solution of the linear 
system. 

Introducing the solution (|19p into equations (|17p we obtain the following equations: 
• the continuity equation, 

I (kxUx + kyUy) + DW = 0, (20) 



• the Navier-Stokes equations. 



Pr 



Ux 



a 

Pr " 



a 

'P'r' 



-W 



the heat equation. 



= -ikxSp + VUx + ikxv'{z)W, 
= -iky5p + VUy + ikyv'{z)W, 
= -D{5p) +VW + v'{z)DW + m, 

ere = {D^ - /c^) e + w, 



(21) 
(22) 
(23) 

(24) 



where the operators Vf = D{vDf ) — vk^ f and Df = df /dz and the wave number k'^ = k'^ + ky 
have been introduced. 

The boundary conditions become 

e=^Ux = Uy = W = Oonz= -1/2; 6 = DUx = DUy = = on z = 1/2. (25) 

Lemma 1 If {Ux,Uy,W,Q,5P), such that {Ux,Uy,W) £ , is a solution of problem S20\\24\ > 
with b.c. ^25\) then it is a solution of the following problem 



Z^v{z) [D^-\k\^]Z +v'{z)-DZ. 



Pr 



Pr 



[D' - \k\^)W = -i?|/.f e - \k\-'rW + D{r{DW)) + \k\^v"{z)W, 



(26) 
(27) 
(28) 



Z = W = DW ^<d = Q, on z ^ -1/2; Z ^ D^W ^0, on z^ 1/2. (29) 

Where Z (z) = Uy (z)ikx — Ux {z)iky, and the velocity functions Ux and Uy can be calculated 
from W and Z. 
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Proof. Pressure is perfectly determined by the velocity {Ux, Uy, W) and the temperature 8, so 
we proceed to eliminate the pressure applying {ik^, iky, D) x {Ux, Uy, W) into the Navier-Stokes 
equations (|21j|23p hence ik^-p^Uy — iky-p;j:Ux leads to (^5]) . 

On the other hand, from the equations (|2m22[) and together the continuity equation (|20p we 
can obtain 

- ^DW ^\k\^SP--r(DW). (30) 
Fr 

Hence one deduces the decoupled equation 

■^Ikl^W = Vlkl^W + Rlkl^e - D \viDW) - -^Dw] - \k\^v"{z)W, (31) 
Fr L Fr J 



(32) 



which rewrites as pS]) and therefore the final system (|26ll28p is obtained. 
Finally, if we define the velocity functions Ux, Uy as, 

Ux = ikx(j} — ikyip, Uy = iky(j) + ikxip, 



using the continuity equation (|20p and knowing that Z (z) = Uy (z) ikx — Ux (z) iky, we obtain 
that (/)(z) and ip{z) are functions of W (z) and Z{z) respectively, 



Ux (z) = I 

Uy (Z) = Z 



, ( DW{z) \ (Z{zy 



(33) 



SO, we can reduce the system (|20ll24p in (|26ll28p without the terms with Ux and Uy. 
3.2 Real eigenvalues 



Theorem 2 The eigenvalues for problem \20\29\) are real. 

Proof. Notice that, for any function / G i7g([— i, i]), we have the identity 



Vffdz ^ ~e I ' v{z) \dhf? + I/I' 



dz 



as well as the identity 



{dl - l)ffdz = - 



M^./l' + |/|' 



dz. 



where dh = hD = D/k = d/dz/k. We consider the equations (|27ll28p and write them as follows, 

cre(z) = P {dl-l)e{z) + W{z), (34) 
-R&{z) = ■^{dl-l)W{z)-k^dl{i^{z)dlW{z))+2k^dh{iy{z)dhW{z)) 

^k^v{z)W{z)-k^dh{y{z)W{z)), (35) 

where we have cleared the temperature function and the derivatives d]^ = ^jD-' = /k^ = 
d^/dz^/k^ . 



J Wedz = a J \Q\^dz + k^ j (|e|^ + |4e|^) dz. 



(36) 
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On the other hand, we consider the complex conjugate in the equation (|35p . we muhiply it by 
W{z) and integrate, obtaining 



/ 



WQdz = j^J(^\W\^ + \dhW\^yz + 



% I Kz) (|W/p + 2\d,W\' + \dlW\' + dim')) dz. (37) 



R 

We thus deduce 

a j lepdz + k' j (lep + |d„e|2)dz = / (iw^p + \dhW\'')dz 



'r 



v{z) {\W\' + 2\dHW\' + \dlW\' + 4(1^^1')) dz. 

(38) 

Considering the imaginary part in relation psp we obtain 

Ima lepdz + / dW^I' + \dHW\'')dz^ = 0. (39) 

So, if (T is not real, we deduce that and W are zero, hence is a trivial solution. The theorem 
is proved. ■ 

3.3 Existence of the eigenvalues 

Before going into the proof of existence of eigenvalues, some previous lemmas are necessary. 
In what follows, we shall make use of the following inequalities: 

Lemma 3 For / e Hl{[-\, \]), we have the Poincare inequality as follows 

- 2 - 

\dkf\'dz>^ r l/pdz, (40) 



v{z)\dhf\'dz>-^ v{z)\f\'dz, (41) 



where Si is a constant depending on ^{z). 

Proof. The inequality (|1(H) is classical. It is obtained through the minimization of \Df\''dz 

^ 2 

under the constraint J^i l/pdz for / element of i?o ([— |i ^]); the solution of this problem being 
given by the function costtz, hence the inequality. On the other hand, we can deduce using (j40p 
and v„iin < v{z) < v„iax the inequality as follows 

Kz)|4/(^)Pdz>^™„ / |d,J(z)pdz > / \f[z)\^dz>^^ I u{z)\f{z)\''dz. 

(42) 

So Lemma |3] is proved. ■ 
Now we introduce a = a/k'^. 

Lemma 4 Vct e J? there exists k E M such that R{a-, k) > 0. 
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Proof. Considering = [—5,5] and using the Lemma 3 and the Holder inequahties we obtain 
from the heat equation (|^7|) 



^+1+ fc2j \MLHn) < \\W\\mn), (43) 
where O = fc^O. From the decoupled equation ^E^, 

+ [-1^%] "("^ + + '^^^"^"^O "^"^o(") - 

(44) 

where R = R/k^. Hence we get the estimate 

(45) 

We must thus deduce from (|45l) that if ct > then 

(l + l) (l + S) llW^lli^(n) <i?l|l^lli^(0), (46) 

hence, using that A; > 



R*(k)={ min iy(z)(k^ + Si) + k^ min (z^(0)) (fc^ + tt^) < i?(/fc), 

V[-l/2,l/2] [-1/2,1/2] ' "J ' 

SO R{k) > because d^(i^(z)) > Vz G [—1/2, 1/2]. On the other hand, if cr < we consider 
from (gg) that 

\a\<min< min i'(z)(l + %], min d?(i^(z)), f 1 + tt 

U-l/2,1/2] V [-1/2,1/2] V 



Define the operator A4 = {1 — dj^ + We have 



and we also obtain that R{k) > 0. 
Lemma 5 Any solution of problem li26Vl29\) satisfies the boundary conditions: 

Proof. First of all, we identify the eigenvalues and eigenvectors of the operator A4 = (l — df^ + a). 
For this purpose, we solve 

(1 + - dl)T = AT 
with boundary conditions r(±^) = 0. This rewrites 

{dl + x-{i + d))e = o. 

As the equation {df^ — uj'^)u = has not non-trivial solution in Hq for lu real, we have necessarily 
A>(1 + (t). 

We introduce uj such that cj^ = A — 1 — ir. We obtain 

T{x) — Acos{LLjkx) + B sm{u!kx), 

from which we deduce the condition to have a non zero non trivial solution sinwfc — 0, hence 
Luk = TiTT. The associated eigenfunction is cos(n7ra:;) for even n and is sin(n7rx) for odd n. We 
denote these functions by Tn(x) and thanks to Fourier analysis, is described by a sum ^ a„T„ 
where X] < +00 and Hq is described by X) a„r„ with ^ n^|a„p < +00. We have M{Tn) = 
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(■^5 — h 1 + (j)T„ = A„((T)r„. Hence the operator A4^^^ is a symmetric operator with the 
eigenvalues ^J,n{^) = \/ ^ni^) and for Q G Hq{[—^, i]) we have A4Q = J2 0'n\/JhiW)Tn- 

As M{a)Q ^WandWe H^, M^^^B = M'^/^W and A^i/^g) e i?o^. Because = 
J2wnTnwithJ2n'^\wn? < +00 and t^nT^ (± f ) = 0. We deduce that A^-i/^iy = Mn ^^^w„T„. 
This is an element of H^, hence we can compute the trace on the boundary, and as the trace on 
the boundary of the cigcnfunctions is zero, we have Al~^/^VF(±i) = 0. Hence M^^^Q{±^) = 
■ 

Theorem 6 Fo?' all k ^ R and for all a > Q, there exists a function R±(k,a) satisfying 
R±{k,a) > R*{k) such that problem I26^28\) with the be i29\) has a unique solution {Z,&,W) € 
(-ffg)^ X Hq, where Hq are the functions of H'^ that fulfill the be for R ~ R±{k,a). This is 
in particular true for a ~ 0. 

Proof. We consider the equations (|27ll28p and we write them as follows: 

7W(5)e(z) = W{z), (47) 
C{d)W{z) = R{k) ■ e(z), (48) 

where the operators Ai and C are defined as 

M = (1-^2+ a), (49) 

C = d,,(Q-^^dH-(^Q-^^+dli,.iz)) (50) 

with Qf = Vf/e = dh (v (z) • dnf) -v{z)- f. 

The equations (|47ll48p lead to following eigenvalue problem 

L{d)Q{z) = R(k)-'^{z) . . 

e(-l/2) = 6(1/2) =0 ^ ' 

where the operator L is not a self-adjoint operator and is a composition expressed as follows 

L = liM. 

Our operator L is not a self-adjoint operator, but £ and M are self-adjoint operators. In fact 
is a positive operator. 



< Twe, e >= (1 + ct) y |e|" + y |d,.er, (52) 

so it is a symmetric operator. As M. is self-adjoint there is a orthogonal basis of i?g with 
the eigenvalues A„ S R. Hence the operator M^^"^ is a symmetric operator with eigenvalues 
yUji = \f\i e R- On the other hand M^l"^ is a closed operator because Sobolev embedding, 
iP' ^ Hq ^ L^. Therefore by Von Neumann theory A4^^^ is a diagonal self-adjoint operator. 
The problem ([51]) is equivalent to the Sturm-Liouville problem as follows 

H{dMz) = R{k)-ip{z) , . 

^(-1/2)^^(1/2)^0 ^^""^ 

where the Dirichlet boundary conditions can be used thanks to lemma [5l the unknown solution 
is If ~ and H = M^^^CM^^^ is a self-adjoint operator because M^^^ and C are self- 

adjoint operators. £ is a self-adjoint operator from Hq to H~^, hence if is a self-adjoint operator 
from Hq to H~^. Therefore R{k) = is not eigenvalue, so the Sturm-Liouville problem 

Hi^Mz) = 

^(-1/2) = ^(1/2) = ^"^^ 
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has only the trivial solution as solution. Using the Frccdholm Alternative theorem we can 
deduce that there exists a function R±{k), which is the smallest eigenvalue of the operator H 
with Dirichlet boundary conditions and the problem (|5ip has a unique and not trivial solution 
&{z). Hence there exists a Green function G{z,£^) such that 



e(z) 

and the theorem is proved. ■ 
3.4 Exchange of stability 



R{k) 



G(z,oe(c)d^, 



Theorem 7 There exists a 3x3 matrix M{R, k) such that the values Ro{k) for which a{Ro{k)^ k) = 
fulfill the condition det(M (Ro{k) , k)) — 0. 

Proof. For any k G Rwe calculate the values of R for which a{R{k), k) = 0. From the equations 
()34ll35p . we consider the heat and decoupled equations as follows, 

I 2 



a 



(55) 
(56) 

(57) 
(58) 



Rk'^e = d{v --^^ DW -k^ (v-^- iy"iz)j W, 
where Vf = D {v (z) • Df) — k^ {v [z) ■ /). Homogenizing equations (|55ll56p we obtain 
{dl - 1) e - ae + = 0, 

RQ = dh (q - dhW -(q-^- 4(i.(z))j W, 

where Qf = Vf/k^, R = R/k^ and 6 = PS. 

Therefore, from the system (|57ll58p . renaming the variables as follows 

Wo = W, (59) 

Wi = dhW^dhWo, 

W2 = u [z) ■ dlW ^ V [z) ■ dhWu 

W^3 = dhW2-2v{z)-Wi, 

6i = dhQ = dhQo, 

we obtain the system of ordinary differential equations for the velocity and the temperature 
functions 



d,,ei = {i + d)-eo-Wo, 



R-eo = dhWs +i^— + „{z) + di{v{z)) ) • Wo 
These equations admit the following matrix formulation 

Wi 
W2 

Oo \Pr 

V ei / 



(z) Pr 



W2 



(60) 
(61) 



dh 






1 











\ 






1 












v{z) 
























2v{z) 





1 








j.{z) + dl{u{z))^ 





(T 





R 





V [z) Pr 

















1 


-1 











1 + CT 


/ 



( Wo\ 

W-i 

W2 
V 01 / 



(62) 
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with boundary conditions 



Wo z 



= Wo[z = 



r 



Qo z 



Qo\z 




(63) 



In order to solve this problem we first consider the initial conditions on z = —1/2, where initial 
conditions for W3 and 81 arc deduced from the equations, 



W2 yz 
i z 



Qi z 



= dhW2 
= dh&o { - 



■ dhWi ( - 

1 

2 
1 
2 



1 



(64) 
(65) 
(66) 



As we do not know the values of dhWi{—l/2), dhW2{—l/2) and (i/i6o(— 1/2) such that the 
boundary conditions on 1/2 are fulfilled, we consider the following general initial conditions 



Wo z 



Wt z = --] = 



W2 z 



2 



-A 



(67) 



1 
2 

W3 ^z = 
Qo (^z 
61 (z = 



and we solve ([62]) with be ^7} for three different sets of values of {A, B, C): (1,0,0), (0,1,0) and 
(0,0,1). We call the respective solutions W\, W2 and W3. The general solution can be written 
as a linear combination of the form 




W (z) = X-Wl{z) + n- W2 {z)+uj- W3 (z) , 



(68) 



where 



Wi {z) = {Wio, Wti,Wi2, Wh, Bio, e^l) (z) , (69) 

We look for a solution of ((62|) with be ([63|l . therefore the solution W must verify the boundary 
conditions on z = 1/2 



(70) 




This system can be expressed as 



M{R, k) 











- 


I) 









(71) 
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where the matrix M{R, k) is defined by 



M{R, k) 




(72) 



Taking = 0, problem (|62j|63p has solution, therefore the linear system ((7T|) has a solution, i.e., 
there is a value of R, Ro{k), such that the system ([70)1 has a solution, i.e., there exist a value 
Ro{k) such that det(M(^o(fc), fc)) = 0. Figure 2 shows this function _F(i?)=det(M(E, fc)), in the 
case of exponential dependence of viscosity on temperature, for 7 = 3- 10^ and a fixed value of 
k, k = 2.54, crossing the axis at Ro{k). ■ 

Corollary 8 For any k £ R there exists values of R for which a > and values of R for which 
a < 0. Therefore, the exchange of stability holds. 

Proof. If we consider ct > 0, problem (|62ll63p has solution, therefore the linear system ((7T|) has 
a solution, i.e., there is a value of R, R'^{k), such that the system (|70p has a solution, i.e., there 
exist a value R~^{k) such that det(M(i?+(/s), k)) ~ 0. And the same applies to ct < 0. ■ 

From here the marginal velocity and temperature fields are easily obtained. For Ro{k), the 
rank of the matrix is less than three and the system to solve is 



A • eio (\) + 1^ ■ e2„ = -uj ■ eso Q 



(73) 



with the free parameter uj. The solution of this system is the growing perturbation after the 
bifurcation. 

For all k Cz R we calculate Ro{k), this function has a minimum, i.e., there is fee such that 
Roikc) = min The real kc is the critical wave number and Rc = Ro{kc) is 

the critical Rayleigh number. As R increases this is the first value of k for which an eigenvalue 
becomes positive. The corresponding mode grows and the conductive state becomes unstable. In 
fact,_jfor any R < Rc, the conductive solution is stable, it looses stability at Rc and it is unstable 
for R > Rc- The curve Ro{k) is the marginal stability curve which is calculated in the next 
section for some values of the parameters. 

4 Numerical Results 

In this section the viscosity is assumed to have an exponential dependence on temperature 

T, 

v{T)^Vo-cM~l{T-T^)), (74) 

where vq is the viscosity at temperature Ti and 7 the exponential rate. Figure 3 shows two 
different viscosity profiles for 7 = 10 and 7 = 3- 10"', that correspond to almost constant 
viscosity and strong variable viscosity, respectively. 

First we have calculated the critical values Ro{k) for each k looking for the first root of 
det(M(i?, k)). The marginal stability curves for 7 = 10 and 7 3 • lO"' can be seen in figure 4. 
The minima in the critical wave numbers kc = 2.68 and Rc = 1082.90 for 7 = 10 and kc = 2.16 
and Rc = 73.51 for 7 = 3- lO'' indicates the critical thresholds for the first bifurcation. In table I 
the critical Rayleigh and wave numbers for several values of 7 are displayed. The corresponding 
critical Rayleigh number, Rc, decreases when 7 increases. Therefore a higher variation of the 
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exponential rate in viscosity favors the instability. The critical wave number also depends on the 
viscosity factor 7. It decreases when 7 increases until 7 = 5- 10"* where it increases again. As 
7 tends to zero, — > 1^0 the results tend to the ones with constant viscosity, kc = 2.68 and 
Rc — 1100.65. Figure 5 shows the corresponding velocity functions of the perturbation field for 
two cases: almost constant (7 ^ 10) and strongly variable viscosity (7 = 3- iC), for different 
Rayleigh numbers, below the critical one, at the critical threshold and above it. In the variable 
viscosity case the velocity is larger in the region where viscosity is smaller, while in the constant 
viscosity case velocity is more distributed along the cell. Results are similar to those in ref. [18] 
where a linear approximation to the exponential dependence is considered. 

5 Conclusions 

We have studied the thermoeonveetive instability problem in an infinite layer under the 
perspective of viscosity as an exponential function of the temperature. We have obtained the 
conductive solution and we have studied its linear stability analysis. We have demonstrated that 
the eigenvalues are real and the existence of eigenvalues and bifurcation. 

With appropriate changes of variable we have considerably simplified the expressions to obtain 
a system of ordinary differential equations which is numerically manageable. We have obtained 
a practical condition which permits easily to calculate the critical bifurcation parameters. 

Finally, we have obtained the stability curves for different values of the exponential rate, 7. 
As the exponential rates decrease the critical Rayleigh number increases. So a higher variation 
of the exponential rate in viscosity favors the instability. There is a non monotone variation of 
the critical wave number with the exponential rate. The vertical velocity for the growing mode 
exhibits more movement in the region where the viscosity is lower. 
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Table captions 



Table I 

Critical wave and Raylcigh numbers {kc,Rc) for different values of 7. 

Figure captions 
Figure 1 
Problem set-up. 

Figure 2 

F {R)=det{M {R, k = 2.54), k ~ 2.54) for 7 = 3 • 10"'; this plot represents a transcritical bifurca- 
tion. 

Figure 3 

Viscosity profiles i'{z) for 7 = 3 • 10'* (dashed hue) and 7 = 10 (solid line). 
Figure 4 

Marginal stability curves for 7 = 3 • 10"* (dashed line) and 7 = 10 (solid line). 
Figure 5 

a) Velocity functions W of the growing mode for different values of the Raylcigh number for 
7 = 10 and k = 2.68; b) velocity functions W of the growing mode for different values of the 
Raylcigh number for 7 = 3 • 10^ and k = 2.16. 
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